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I .  INTRODUCTION 


The  simplest  theory  for  nonlinear  wave  propagation  in  a  straight, 
uniform  rod  may  be  written  as  follows. 


S 


Z 


pv. 


V 


Z 


(1) 


The  axial  stress,  S,  depends  only  on  the  axial  strain,  e,  the  axial 
particle  velocity  is  v,  and  the  reference  density,  p,  is  a  constant.  The 
axial  coordinate  along  the  rod  is  Z,  t  is  time,  and  subscripts  denote 
partial  differentiation.  In  this  version  of  rod  theory,  which  is 
formally  the  same  as  one  dimensional,  inviscid  gas  dynamics,  the  char¬ 
acteristic  speed  of  propagation,  c,  is  given  by  pc2  *  dS/de,  and  dis¬ 
continuities  in  stress,  strain,  or  velocity  . all  travel  at  the  shock 
speed  given  by  pc2  =  [S]/[e],  where  [.]  denotes  the  jump  in  a  quantity 
across  the  shock. 


Equation  (1)  has  been  used  to  interpret  wave  propagation  experiments 
in  which  the  strain-time  profile  is  measured  at  several  stations  in  the 
r.od.l»2  If  each  level  of  strain  propagates  at  a  constant  velocity, 
c  (e) ,  then  both  stress  and  particle  velocity  may  be  computed  as  a  func¬ 
tion  of  strain  as  follows. 


re(Z,t) 

v  =  v  -  I  c  de 

0  Je  P 

o 

(2) 

S  -  s0  +  J e  PCp  de 


The  initial  values  ahead  of  the  wave  are  vQ,  and  eQ.  Equation  (2)  is 
a  simple  wave  solution  to  (1) ,  exhibiting  amplitude  dispersion,  but  not 
geometric  dispersion. 


\j.  F.  Bell,  The  Experimental  Foundations  of  Solid  Mechanics,  Handbuch 
der  Physik,  Vol.  VIa/1,  Springer- Verlag,  New  York,  1973. 

2 

G.  E.  Hauver,  "Penetration  with  Instrumented  Rods,"  Int.  J.  Eng.  Sci, 
16,  1978,  pp.  871-877. 
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A  rod,  no  matter  how  slender,  is  actually  a  three-dimensional 
object,  of  course.  The  dispersive  influence  of  a  finite  diameter  is  well 
known  in  linear  theories,  eg. 3*4  it  is  less  well  known  in  nonlinear 
theories5*6  where  it  may  be  expected  to  play  a  role  as  well,  whenever 
the  length  scale  in  a  wave  pulse  is  comparable  in  magnitude  to  the  rod 
diameter.  In  this  paper  the  effects  of  finite  lateral  dimensions  will 
be  considered  by  modeling  a  rod  as  a  one-dimensional  elastic  structure 
with  one  internal  variable  used  to  represent  the  transverse,  axisymmetric 
motion.  Such  a  structure  is  a  special  case  of  an  intrinsic  rod  theory 
as  described  by  Antman^. 


II.  ONE-DIMENSIONAL  EQUATIONS  AND  RELATIONSHIP  TO 
THREE-DIMENSIONAL  ELASTICITY 

A  straight  cylindrical  rod  of  radius  a  is  assumed  to  have  an  elastic 
stored  energy  density  per  unit  length 

W  -  ir  a2  W(w',u,u').  (3) 

The  axial  displacement  is  w,  u  is  a  measure  of  radial  strain,  and  the 
dash  represents  differentiation  with  respect  to  the  axial  coordinate,  Z. 
The  kinetic  energy  density  per  unit  length  is  given  by 

K  ■  tt  a2^sPjW2  +  %p 2u2) ,  (4) 

where  and  P2  are  the  appropriate  mass  densities  (made  definite  by 
equation  (13))  for  the  axial  and  radial  motions  respectively,  and  the  dot 
denotes  differentiation  with  respect  to  time.  The  Euler-LaGrange  equa¬ 
tions  corresponding  to  (3)  and  (4)  are  given  by  (5) . 


T 

R.  Skalak,  "Longitudinal  Impact  of  a  Semi-Infinite  Circular  Elastic  Bar," 
J.  Appl.  Mech.,  24,  1957,  pp.  59-64. 

^W.  A.  Green,  "Dispersion  Relations  for  Elastic  Waves  in  Bars,"  in 
Progress  in  Solid  Mechanics,  Vol.  I,  ed.  I.  N.  Sneddon  and  R.  Hill, 

North -Ho 11 and,  Amsterdam,  1960. 

5J.  H.  Shea,  "Propagation  of  Plastic  Strain  Pulses  in  Cylindrical  Lead 
Bars,"  J.  Appl.  Phys.,  39,  1968,  pp.  4004-4011. 

6G.  P.  DeVault,  "The  Effect  of  Lateral  Inertia  on  the  Propagation  of 
Plastic  Strain  in  a  Cylindrical  Rod,"  J.  Mech.  Phys.  Sol.,  13,  1965, 
pp.  55-68. 

7S.  S.  Aritman,  The  Theory  of  Rods,  Hanbuch  der  Physik,  Vol.  VIa/2, 
Springer-Verlag,  New  York,  1972. 
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(5) 


S'  -  Pjtt 
q'  -  P  »  p2u 

The  forces  S,  P,  and  Q  are  obtained  from  the  stored  energy  as  follows. 


S 


3W  n  3W 
3w'  ’  P  “  3u 


(6) 


Equations  (5)  and  (6)  may  be  interpreted  in  terms  of  the  three- 
dimensional  theory  of  nonlinear  elasticity.  For  the  axisymmetric  motions 
considered  here,  the  deformation  from  an  unstressed  reference  configura¬ 
tion  (R,0,Z)  to  the  present  configuration  (r,0,z)  may  be  written 


r  =  r(R,Z,t),  6  =  0,  z  e  z(R,Z,t)  ,  (7) 

and  the  corresponding  equations  of  motion  are  given  by 


aTrZ  TrR  _  Te© 
3Z  R 


pr 


3T 


zR 


8T 


zZ 


nzR 


3R 


3Z 


pz 


(8) 


The  stress  components  used  here  are  obtained  from  the  unsymmetric  Piola- 
Kirchhoff  tensor  referred  to  unit  basis  vectors  in  the  given  cylindrical 
coordinates  (r,0,z)  and  (R,0,Z). 

T  *  Tx“  9i  a  eo  (9) 

If  (8)..  is  multiplied  by  R,  and  if  the  resulting  equation  and  (8)2  are 
both  averaged  over  the  cross-section,  equations  (10)  are  obtained  for  the 
case  of  stress-free  lateral  surfaces. 


i 

4 


1 


-4  si  p11 « 

ira  * 

dA  -  (jTR 


+  T00)dA 


Now  if  the  following  identifications  are  made 


v 


7 


S  ~  fllZ  dA  , 
ira^ 


p 

1 

7?-> 

f  (TrR  ♦  T90)  dA  , 

(11) 

Q 

-T  J 

iTa 

f RTrZ  dA  , 

and  if  a  first  approximation  for  the  position  functions  r  and  z  is 
assumed  to  be 


r  -  R[1  +  u(Z,t) ]  , 
z  =  Z  +  w(Z,t)  , 


(12) 


then  equations  (5)  are  obtained  with 

2 

pl  =  p  »  p2  = 


(13) 


In  all  of  the  preceding  discussion  it  has  been  tacitly  assumed  that 
the  Z-axis  is  an  axis  of  material  symmetry  so  that  no  angular  motion  will 
occur  anywhere  in  the  cross-section.  Furthermore,  since  the  strain 
energy  must  be  invariant  under  reversal  of  the  z-axis,  and  if  it  is 
assumed  to  be  invariant  under  reversal  of  the  Z-axis  as  well,  W  must  have 
the  property 


W(w',u,u')  =  W(w',u,-u') 


(14) 


It  follows  that  Q  is  odd  in  u'  and  S  and  P  are  even  in  u' . 

It  should  be  remarked  that  equations  (5)  and  (6)  have  been  used  to 
describe  one-dimensional  waves  in  porous  materials°»9>10,  and  it  has 


O 

J.  W.  Nunziato  and  E.  K.  Walsh,  "On  the  Influence  of  Void  Compaction  and 
Material  Non-uniformity  on  the  Propagation  of  One-Dimensional  Accelera¬ 
tion  Waves  in  Granular  Materials,"  Arch.  Rat.  Mech.  Anal.,  64,  1977, 
pp.  299-316  and  Adendum,  Arch.  Rat.  Mech.  Anal.,  67,  1977,  pp.  395-397. 

Q 

J.  W.  Nunziato  and  E.  K.  Walsh,  "One-Dimensional  Shock  Waves  in  Uni¬ 
formly  Distributed  Granular  Materials,"  Int.  J.  Solids  and  Structures, 
14,  1978,  pp.  681-689. 

10D.  F.  Parker  and  B.  R.  Seymour,  "Finite  Amplitude  One-Dimensional 
Pulses  in  an  Inhomogeneous  Granular  Material,"  Arch.  Rat.  Mech.  Anal., 
72,  1980,  pp.  265-284. 
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also  been  suggested  that  perhaps  they  could  be  used  to  describe  waves  in 
a  layered  composite  .  There  are  undoubtedly  other  applications  as  well. 


III.  STATIC  SOLUTIONS 

It  is  convenient  to  distinguish  two  types  of  static  solutions.  In 
the  first  type  u'  is  identically  zero,  so  equations  (5)  reduce  to 


S(w',u,0)  ■  SQ  , 
P(w',u,0)  =  0  . 


(15) 


Since  the  reference  configuration  has  been  assumed  to  be  unstressed  >  it 
is  required  that  when  w*  «  u  «  0,  (15)  is  satisfied  with  S  =  0.  It  will 
be  further  assumed  that  the  slope  along  the  curve  defined  by  (15)  2  is 
always  negative  and  bounded. 

.  P  f 

-  >  -  -  jr-  =  v(w',u,0)>  °  (16) 

u 

where  the  subscripts  u  and  w*  denote  partial  differentiation  with  respect 
to  the  arguments  of  P.  This  corresponds  to  the  physically  reasonable 
assumption  that  the  rod  contracts  (expands)  laterally  as  it  is  stretched 
(compressed).  Finally,  it  will  be  assumed  that  S  increases  monotonically 
as  w*  increases  along  the  curve  P(w*,u,0)  =  0.  That  is 

Sw,  -  vSu  =  E(w»,u,0)  >  0.  (17) 


The  bar  modulus,  E,  and  Poisson* s  ratio,  v,  will  be  considered  as  defined 

by  (17)  and  (16)  with  the  inequalities  holding  as  well  for  all  values  of 

their  arguments,  not  just  along  the  curve  of  (15) 2 *  Since  W  has  a  minimum 

at  the  unstressed  reference  configuration  (w*  =  u  =  u*  =  0) ,  Pu  >  0 

there,  so  by  (16)  P  f  >  0  and  P  >  0  everywhere.  This  also  implies  from 
w  u 

(17)  that  W  t  fW  -  W  ,  >  0  so  that  necking  instabilities  have  been 
v  *  w*w*  uu  u  w*  ® 

ruled  out  (see  Antman  [12],  pp.  97). 

In  the  second  type  of  static  solution  u*  is  not  identically  zero. 
There  are  two  integrals  of  (5) . 


F.  McCarthy,  private  communication. 

12 

S.  S.  Antman,  "Qualitative  Theory  of  the  Ordinary  Differential  Equa¬ 
tions  of  Nonlinear  Elasticity,"  in  Mechanics  Today,  V.  1,  ed.  S.  Nemat- 
Nasser,  Pergamon,  New  York,  1972. 
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(18) 


S(w',u,u')  -  SQ 

w'QCw'.u.u’)  +  w’So-W(w' ,u,u')  »  B 

If  these  may  be  solved  for  u'  and  w'  as  functions  of  u  and  the  constants 
S9  and  B,  then  u  may  be  found  by  quadrature,  and  the  axial  strain  w'  is 
given  parametrically  through  u.  The  study  of  equations  (18)  will  not  be 
pursued  further  here  except  to  remark  that  the  version  arising  from 
linear  elasticity  gives  boundary  layer  solutions  with  Q  decaying  exponen¬ 
tially  from  the  ends.  Phase  plane  analysis  indicates  that  the  behavior 
of  the  fully  nonlinear  equations  is  similar. 


IV.  SOLUTIONS  OF  THE  LINEARIZED  EQUATIONS 

Before  considering  the  dynamic  nonlinear  equations,  it  will  be  use¬ 
ful  to  examine  some  of  the  properties  of  solutions  of  the  linearized 
equations.  The  three-dimensional  strain  energy  for  an  isotropic  mater¬ 
ial  with  displacement  field  given  by  (12),  when  averaged  over  the  cross- 
section,  leads  to  (19)  for  the  forces. 


S  =  (A+2y)w'  +  2Au  , 

P  =  4(A+p)u  +  2Aw’  ,  (19) 

r  ,2  1 

Q  •  ha  vu  > 


where  A  and  p  are  the  usual  Lame  constants.  With  nondimensional  variables 
C  =  Z/a  and  t  =  ct/a,  where  c  is  a  speed  that  is  unspecified  for  now,  and 
with  w  =  w/a,  equations  (5)  become 


w«  +  2  x^T 


u 


u  -  1^8  u  +  4  —  w  1 
«  l  v  VI  C  J 


—x-  W 

.  2  tt 
'1 

2 

V  u 
,2  TT 


(20) 


where  subscripts  again  denote  partial  differentiation.  Subscripts 
indicate  partial  differentiation,  and  the  overbar  has  been  dropped  from 
w.  The  longitudinal  and  shear  wave  speeds  are  c^  ■  /(A+2y)/p  and 

c,  *  *Wp  .  These  equations  in  dimensional  form  were  first  given  by 
Mlndlin  and  Herrmann13. 


13R.  D.  Mindlin  and  G.  Herrmann,  "A  One-Dimensional  Theory  of  Compres- 
sional  Waves  in  an  Elastic  Rod,"  Proceedings  of  the  First  U.S. 
National  Congress  of  Applied  Mechanics,  1950,  pp.  187-191. 
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L 


Since  all  the  coefficients  are  constants,  either  w  or  u  may  be 
eliminated  from  (20)  to  obtain  (21). 


2  .  2 
c^  3t 


c  a 

2  ^  2 
3t 


X+y  Cb  ( 3 2  c2  3 2 

y  2  L  2  2  2 

C1  %  3t 


(w,u)  *  0 


The  one-dimensional  bar  speed  is  c,  =  /E/p~  where  E  is  Young's  modulus, 

E  =  y(3A  +  2u)/(A+p).  Equation  (2 1J  exhibits  what  Whitham  has  called  a 
hierarchy  of  wave  speeds^*  15.  Whereas  discontinuities  can  only  propagate 
with  the  speeds  c^/c  or  C2/C,  it  is  well  known  that  in  a  thin  bar  the 
main  disturbance  in  a  pulse  travels  with  speed  c^/c.  Whitham  has 
discussed  extensively  the  hierarchical  case  where  the  orders  of  the 
highest  and  lowest  derivatives  differ  by  one.  Wu*°  has  pointed  out  that 
the  case  where  the  orders  differ  by  two,  as  here,  also  exhibits  a  stable 
hierarchical  wave  structure,  so  long  as  c2<cb<cly  w^ich  is  always  the  case 
for  linear  elastic  rods. 

In  Whitham's  case  the  higher  order  waves  decay  exponentially,  and 
the  lower  order  wave,  which  carries  the  main  pulse,  has  a  diffusing  front. 
In  the  present  case,  discontinuities  in  the  higher  order  waves  do  not 
decay  with  time,  but  the  pulse  decays  rapidly  behind  the  front.  The  main 
pulse  travels  with  the  lower  order  speed,  and  it  too  has  a  diffusing 
front . 

To  give  substance  to  the  preceding  remarks,  consider  a  boundary  - 
initial  value  problem  for  (20) . 

B. V.  :  u  (0,t)  =  0,  w  (0,t)  =  f(r) 

Q  (22) 

I.V.  s  u(C,0)  =  w(C,0)  =  0 

A  straightforward  application  of  the  LaPlace  transform  on  t  leads  to  the 
following  representations  for  u  and  w. 


G.  B.  Whitham,  "Some  Comments  on  Wave  Propagation  and  Shock  Wave  Struc¬ 
ture  with  Applications  to  Magnetohydrodynamics Comm,  on  Pure  and 
Appl.  Math,  12,  1959,  pp.  113-158. 

^G.  B.  Whitham,  Linear  and  Nonlinear  Waves,  John  Wiley,  New  York,  1974. 


T.  T.  Wu,  "A  Note  on  the  Stability  Condition  for  Certain  Wave  Propaga¬ 
tion  Problems,"  Comm,  on  Pure  and  Appl.  Math.,  14,  1961,  pp.  745-747. 
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?(p)  is  the  transform  of  f(x),  and  m^(p)  andm2(p)  are  roots 
negative  real  parts  (for  outgoing  waves)  of 


2  2c 

m  -P  — 

v  c2 


2  2 
m  -p 


(25) 


In  each  equation  of  (23)  the  first  integral  represents  a  wave  whose  front 
travels  at  speed  c^  and  the  second  integral  represents  a  wave  whose  front 
travels  at  speed  c*. 

•  To  examine  the  behavior  of  the  solution  near  the  wavefronts,  let  p 
be  large.  Expansions  for  m^  and  m2  are 
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m 


1 


m 


2 


(26) 


c 


The  leading  terms  in  the  integrals  for  w  may  be  written  asymptotically 
as 


In  (27)  the  square  root  terms  are  equivalent  to  or  m2  within  0 
when  r^2  and  r22  are  given  by 


2  2 

r2  fl"1  _Cb 

rl  8  - 2 


2  2 
2  .  *b  -  C2 

r2  =  8  - 2 - 


(28) 


The  first  integral  in  (27)  vanishes  for  x  <  —  c  and  the  second  vanishes 

C1 

for  x  <  —  c-  Thus,  near  the  leading  wavefront  only  the  first  integral 
c2 

is  required,  which  may  be  written,  after  a  little  manipulation  and  use 
of  the  convolution  theorem,  as  follows. 
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JQ  and  are  Bessel  functions, 
is  the  Heavisde  step  function  and 


variable  y  = 
ence  17,  pg.  64) 


-  nf 2  -  4 


Suppose  that  f(x)  =  w  0  h(x)  where  h(x) 
i  w _  is  a  constant.  The  change  of 


C  and  use  of  the  identity  (taken  from  Refer- 


jx(y)  dY 
J  2  "  2 

Vy  +  n 


(30) 


now  reduces  (29)  to 


w 

x 


w 


o 


X 


+  r 


JjCy)  dY 


•  (31) 


Equation  (31)  shows  that  at  points  removed  from  the  boundary  the  velocity 
becomes  exponentially  small  shortly  after  passage  of  the  wave  front. 
Equation  (29)  is  plotted  in  Fig.  1  with  f(x)  =  h(x). 

Similar  analysis  may  be  applied  to  the  other  integrals  in  (23). 

The  second  integral  for  w^  represents  a  wave  whose  front  moves  with 
speed  c2,  and  is  similar  in  form  to  (29),  but  smoothed  by  a  double 
integration  in  x.  u  has  a  contribution  at  each  wavefront  that  is 
similar  to  (29),  but  Smoothed  by  a  single  integration  in  x. 


^F.  Bowman,  "Introduction  to  Bessel  Functions,"  Dover,  New  York,  1958. 
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To  examine  the  main  pulse  in  (23)  it  is  convenient  to  follow 
Whitham's  approach  to  the  problem  (Reference  15,  pg.  345).  First,  the 

exponent  in  one  of  the  integrals  is  written  as  t^p  +  m ^ (p)  -p)  ,  and  then 

the  saddle  point  method  is  used  to  find  an  asymptotic  expansion  valid 

at  large  t.  For  a  fixed  value  of  C/x  saddles  occur  at 

% 

3p  (p  ♦  VP)  y)  •  1  ♦  V(p)  T =  °*  c32) 

Solution  of  (32)  gives  p  =  p  (-p)  »  the  exponent  is  maximized  with 
respect  to  C/t  when  m^[p(C/T)]  =  0  because  of  (32)  and  (33). 


diW  [p  +  micft  t]  =  mi(p) +  f1  +  v  (p>  f]  d v!hr =  0  (33) 

2 

From  (25)  if  =  0,  either  p  =  0  or  p  *  -8  — ^  — j-  •  By  differentiating 

c 

(25)  once  with  respect  to  p,  it  is  clear  that  the  correct  choice  is  p  =  0, 
and  by  differentiating  (25)  a  second  time  with  respect  to  p,  it  is  found 
that 

mj'CO)  =  ±  §•  ,  (34) 


and  so  from  (32) 


(35) 


Equation  (35)  shows  that  the  main  disturbance  travels  with  the  bar 
speed,  c^.  The  correct  asymptotic  expansions  will  be  obtained  by  expand 
ing  the  integrands  in  (23)  about  p  =  0.  The  following  discussion  will 
be  simplified  by  taking  c  =  c^. 

It  has  been  shown  that  m1(0)  =  0,  but  it  does  not  follow  that 

2~  2  4 

n^fO)  =  0.  Thereforg,  let  m  =  a  +  bp  +  dp  +  ...  in  (25).  By 
equating  powers  of  p  it  is  found  that 


(cl2  '  °b2)  (°b2  -  « 2 *) 


(36) 
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Choosing  roots  to  correspond  to  outgoing  waves  or  bounded  functions  at 
infinity  gives 


m 


1 


3 


P 


+ 


> 


m 


2 


(37) 


Since  the  second  integrals  for  u  and  w  in  (23)  give  no  contribution 

c 

until  t  >  —  c,  near  the  head  of  the  main  pulse  the  velocity  is  given 
2 

asymptotically  as  follows 


w 

T 


1 

2iri 


*^r 


P(t-;)  - 


;(i  *  0(p2)) 


dp 


where  d^ 


l 

8 


(ci2  -  %2)  {%2  - 


(38) 


By  the  convolution  theorem,  (38)  may  be  written 


where  Ai(-)  is  the  Airy  function  with  integral  representation  given  in 
Abromowitz  and  Stegun*  ,  formula  10.4.32.  For  a  step  function  input, 
f(-r)  =  h  (t),  w  is  simply  the  integral  of  the  Airy  function  and  is 
exponentially  small  for  £  >  r  and  oscillates  about  an  amplitude  of  1  for 
t  >  £.  Fig.  2  shows  the  waveform  for  a  step  input.  Fig.  3  shows  the 
composite  waveform  schematically  for  both  the  leading  wave  and  the  main 
pulse.  Note  the  progressive  separation  and  decreasing  importance  of  the 
leading  fast  wave. 
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M.  Abromowitz  and  I.  A.  Stegun,  Handbook  of  Mathematical  Functions, 
Dover,  New  York,  1965. 
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Figure  2:  Dispersion  in  Main  Pulse.  Linear  Case. 


Composite  Sketch  Showing  Separation  of  Leading  Wave  and  Main  Pulse.  Linear  Case 


The  initial  development  of  this  separation  was  shown  numerically  by 
Nunziato  and  Walsh  in  their  consideration  of  small  amplitude  waves  in 
porous  soils***  where  the  equations  are  formally  the  same  as  (20). 


V.  NONLINEAR  WAVE  EQUATIONS 


The  multiple  wave  structure,  which  the  linearized  equations  show  so 
clearly,  has  an  exact  counterpart  in  the  fully  nonlinear  theory.  The 
characteristic  wave  speeds  of  (5),  which  carry  discontinuities  in  w' ' 
and  w  or  u"  and  ii ,  are  given  by 


=  Swt(w',u,u') 
p2c22  =  (w’ *u*u') 


(40) 


On  the  other  hand,  the  simplest  one -dimensional  theory  for  a  rod,  which 
is  given  by  (1) ,  and  which  may  be  obtained  from  (5)  by  setting  u'  =  0, 
P(w',u,0)  =  0,  has  characteristic  speeds  given  by 

■"S2  ■  E  <  v  .  <41> 

where  E  is  defined  in  (17) .  In  analogy  with  the  linear  theory  it  might 
be  expected  that  (41)  represents  a  low  order  approximation  for  wave 
speeds  in  the  main  pulse,  but  that  dispersive  effects  due  to  the  finite 
diameter  will  modify  the  pulse  shape.  Since  the  full  rod  theory,  (5), 
is  nonlinear,  it  is  to  be  expected  that  wave  speeds  will  be  modified  as 
well. 


Although  transform  methods  cannot  be  used  on  the  full  equations,  it 
is  possible,  by  suitable  scaling  and  choice  of  independent  variables,  to 
find  asymptotic  results  for  the  various  wave  orders  directly  from  the 
equations  themselves  rather  than  from  representations  for  exact  solutions. 
The  method  of  relatively  undistorted  waves,  introduced  by  Varley  and 
Cumberbatch^®  has  been  used  recently  by  Seymour  and  Parker*®  to  examine 
the  leading  fast  wave  in  (5) .  The  same  approach,  also  called  the  method 
of  modulated  simple  waves,  may  be  used  to  investigate  the  structure  of 
the  main  pulse. 

To  simplify  the  analysis,  it  will  be  assumed  that  the  stored  energy 
in  (3)  has  the  decomposition 


E.  Varley  and  E.  Cumberbatch,  "Non-Linear  High  Frequency  Sound  Waves," 
J.  Inst.  Maths.  Applies.,  2,  1966,  pp.  133-143. 
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W(w' ,u,u') 


(42) 


-  WjCw'.u)  ♦  W2(u') 


Let  the  dependent  variables  be  u,  v  *  w,  and  e 
are  rewritten  as 


Equations  (5) 


S'  -  P*  , 
0*  -  P  m  P  , 
v'  *  e  . 


(43) 


First  define  a  new  set  of  independent  variables 


C  =  6  Z,  x 


t. 


(44) 


where  6  is  a  small  parameter.  This  scaling  implies  that  fixed  values  of 
C,  t  will  correspond  to  long  times  and  large  distances  as  6  tends  to  zero. 
Next  define  a  second  set  of  independent  variables  in  terms  of  c ,  t . 


a  =  o(C,t) 
z  =  6? 


(45) 


In  terms  of  these  new  variables,  derivatives  with  respect  to  Z,  t  are 
replaced  as  follows. 


J_  _  A  „  _A  +  *3/2  _3_ 
3Z  "  6  C  3a  6  3z 


(46) 


_3 

3t 


**S  3 

6  a  — 
T  3a 


In  (45)  a  is  a  "fast"  and  z  a  "slow"  variable.  The  reason  for  the  odd 
combination  of  scaling  in  (44)  and  (45)  will  become  apparent  later.  It 
turns  out  that  a  has  properties  like  a  characteristic  variable.  Accord¬ 
ingly,  let  b>  and  k  be  defined  as 


“t  *  « 


(47) 


a  is  constant  along  curves  with  slope 


dc 

3? 


»e 


(48) 


It  is  convenient,  instead  of  u  and  tc  to  use  the  slowness,  s  ■  —  ,  and 

1  c 

the  incremental  arrival  time  l  -  —  .  In  terms  of  these  variables,  the 
compatibility  condition  becomes 
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(49) 


6  *2  . 


and,  taking  (42)  into  account,  equation  (43)  becomes 

sSuu0  *  plv«  +  *Sw ' eo  *  **{Suuz  +  Sw'Ez} 


P  u  -  sP  ,v 
u  o  w'  a 


*  6iJ- P  ,v  +  0  ,u  _  -  p-u  \ 
|  W'  Z  Tl'  T £?  2  TTTr 


(50) 


sv  +  c  *  6l\ 
a  a  z 


where  derivatives  in  t  and  t  are  retained  on  the  right  hand  side  of  (50)2 
for  compactness. 


To  lowest  order  the  right  hand  sides  of  (50)  may  be  set  equal  to 
zero.  This  requires  that  the  determinant  of  coefficients  of  ua,  v  ,  e 
must  vanish  giving  a 


1 


(51) 


and  u  ,  v  ,  e  must  occur  in  the  ratio 
a  a  a 

u  :  v  :e  =  -  v:  -  c,:  1  .  (52) 

a  a  a  d 

With  the  energy  partition  of  (42)  neither  S  nor  P  depends  on  uf.  Thus 
(52)  may  be  integrated  to  find  two  of  u,v,e  in  terms  of  the  third.  For 
example 


i 


V  =  v(u) 
e  -  e(u) 


(S3) 


The  right  hand  side  of  (50)  must  satisfy  a  compatibility  condition, 
multiplication  of  (50)  by  the  left  proper  vector 


Thus 


(-1,  sv,  sSw,) 


(54) 


gives  the  transport  equation 


S  -  sEv 
z  z 


sv  p. 


r-l 

CM 

3 

rs  3/s 

2  * 

3a 

Ll'Sa  (l  Ua)J 

-  i  d  Gr  sKi  ua)]j* 


(55) 


left  hand  side 


"•  »  0  ,.  Since  S,s,E  and  v  all  depend  on  u  through  (53),  the 
side  of  (55)  may  be  written  as  $(u)uz  so  that  the  dependent 
variable  is  u  alone.  However,  the  mapping  a(c>t)  must  be  developed 
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simultaneously  with  the  solution  u(a,z)  since  the  incremental  arrival 
time  i  is  also  unknown. 


The  properties  of  (55)  will  not  be  studied  further  here  except  for 
the  following  two  remarks.  First,  the  linearized  version  of  (55)  can  be 
shown  to  have  solutions  like  (39),  and  second,  for  small  (but  finite) 
amplitude  waves,  the  transport  equation  reduces  to  the  Korteweg-de  Vries 
equation.  This  latter  fact  may  also  be  deduced  directly  from  (5)  follow¬ 
ing  the  prescription  of  Leibovich  and  Seebass  .  Many  other  authors  have 
also  noted  that  the  Korteweg-de  Vries  equation  describes  longitudinal 
waves  in  thin  rods,  eg.,  Nariboli  and  Sedov.  2 


VI.  STEADY  WAVES 


Equation  (49)  indicates  that  if  i,  =  0(1),  then  s,  and  hence  c^,  is 
almost  constant.  Therefore,  some  solutions  of  (55)  should  be  nearly 
steady  waves,  perhaps  approaching  them  asymptotically  in  time.  Accord¬ 
ingly,  the  properties  of  steady  solutions  to  (5)  will  be  examined. 

A  steady  wave  is  one  in  which  the  solution  to  (5)  may  be  written 


w  =  w(Z  -  ct) 
u  •  u(Z  -  ct) 


(56) 


where  c  is  an  unspecified,  but  constant,  wave  speed, 
become 


S' 


w" 


Q'  -  P 


» 


The  equations 


(57) 


where  the  prime  now  denotes  differentiation  with  respect  to  Z  -  ct. 


2 1 

S.  Leibovich  and  A.  R.  Seebass,  "Examples  of  Dissipative  and  Dispersive 
Systems  Leading  to  the  Burgers  and  the  Korteweg-de  Vries  Equations," 
in  Nonlinear  Waves,  S.  Leibovich  and  A.  R.  Seebass,  eds.,  Cornell 
University  Press,  Ithaca,  N.Y.,  1974. 
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G.  A.  Nariboli  and  A.  Sedov,  "Burgers ’s-Korteweg-de  Vries  Equation  for 
Visscoelastic  Rods  and  Plates,’1  J.  Math.  Anal*  and  Appl,,  32,  1970, 
pp.  661-667, 


There  are  two  integrals  of  the  motion. 


S  ■  PjC  w'  +  A 


Q  u'  ♦  A  w*  -  W 

■  I '2(p2“’2  -  *  8 


(58) 


where  A  and  B  are  constants.  Thd  second  integral  is  found  by  multiplying 

(57).  by  w'  and  (52)?  by  u',  adding  the  two  equations,  and  noting  that 

A  “  /  3W  \  1  12/  2  2\  * 

the  result  may  be  written  as  u'  +  g^jr  w'  -  W)  =  ^  c  (pjW  -t^u'  /  . 

After  integration  (58) ^  is  substituted  into  this  expression  to  give  (58)2 


If  (58)  can  be  solved  for  w'  and  u' ,  they  may  be  expressed  as 
follows. 


w'  =  F (u;  A,B,c) 
u'2  =  G(u;  A,B,c) 


(59) 


2 

Equation  (59)2  is  written  in  terms  of  u'  rather  than  u'  since  both  of 

(58)  are  even  in  u' .  Now  (59) ^  may  be  solved  by  quadrature  and  then 

(59) .  gives  the  distribution  in  longitudinal  strain  parametrically 
through  u. 


To  examine  the  qualitative  nature  of  solutions  and  the  role  of  the 
constants  A,  B,  and  c,  first  consider  the  case  A  =  B  =  0.  For  ease  of 
analysis  the  strain  energy  will  be  assumed  to  be  partitioned  as  in  (42) . 
For  small  but  finite  amplitude  waves  let  N  be  expressed  as  a  power  series 
First  note  that 


where  W  is  the  three-dimensional  strain  energy  per  unit  volume.  With 
(12)  used  in  (60) ,  W  may  be  expressed  as 

W  =  j  (A+2y)  w'2  +  2\uuw'  +  2(x+u)u2 

+  j  a2  uu'2  +  C^u3  +  C2U2w'  (61) 

+  C.uw'2  +  C.w'3  +  _  , 

3  4 
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where  the  C's  are  combinations  of  higher  order  elastic  moduli. 


Equation  (58) ^  with  A  =  0  gives 


w'  = - — - j  u  +  °(u2)  > 

X+2y  -  pc 

and  (58) 2  with  A  =  B  =  0  and  the  use  of  (62)  gives 

2  2  3 

uf  =  Cu  +  Du  +  ...  , 

where 

r  -  CE-pc2) 

^  ~  2  2  2  ’ 
a  (y-pc  ) (X+2y-pc  ) 


(62) 


(63) 


(64) 


and  D  may  be  either  positive  or  negative,  depending  on  the  exact  combina¬ 
tion  of  higher  order  elastic  moduli.  C  may  be  either  positive  or  negative 
as  well. 


C  >  0  if  c  <  c  <  c. 

b  1 

C  <  0  if  c2  <  c  <  cfe 


(65) 


where  c^,  c^,  and  are  the  linear  elastic  wave  speeds  in  this  context. 

2 

Figure  4  shows  sketches  of  uf  vs.  u  and  u?  vs.  u  for  the  various 
combinations  of  signs  of  C  and  D.  The  closed  loops  in  (4e)  and  (4f) 
correspond  to  solitary  waves  that  are  symmetrical  in  Z  about  the  point  of 
maximum  amplitude  and  extend  from  Z  -  -«>  to  Z  *  +»  .  The  open  curves  in 
(4e)  -  (4h)  correspond  to  waves  in  which  the  amplitude  grows  without 
limit.  Clearly,  to  be  useful,  solutions  of  this  kind  must  join  onto 
unsteady  solutions  before  the  amplitudes  become  too  extreme. 


Now  consider  the  case  where  A  =  0  but  B  t  0.  Instead  of  (63)  the 
following  holds. 

u* 2  =  B  +  Cu2  +  Du^  (66) 

The  constant  B  is  proportional  to  3,  but  has  opposite  sign.  B  merely 
shifts  the  curves  in  (4a) -(4d)  vertically.  Figure  5  shows  two  of  the 
possibilities.  The  closed  loops  in  (5c)  and  (5d)  correspond  to  periodic 
solutions. 


Next  consider  the  case  when  both  A  i  0  and  B  /  0.  The  integrals  of 
the  motion  may  be  rewritten  as  follows 


s  -  So  «  Plc2(w'-w;) 

Qu'  ♦  SQ(w'-w^)  -  (W  -  WQ)  (67) 

'I'2!' 2U'2  '  8  • 

where  WQ  »  W(w',u  ,0),  etc.  These  integrals  are  the  same  as  before,  but 
now  they  are  centered  on  the  point  (w^,uo,0)  and  SQ  plays  the  role  of  A. 
The  character  of  the  solutions  is  exactly  the  same  as  in  Figures  4  and  5, 
but  now  the  curves  are  shifted  to  the  right  (SQ<0)  or  left  (SQ>0) . 

Finally,  suppose  that  the  elastic  constants  are  such  that  D  =  0. 
Instead  of  (66) ,  equation  (68)  must  be  used. 

u'2  =  B  +  Cu2  +  Hu4  (68) 


This  should  occur  when  the  expansion  is  centered  at  an  inflexion  point  of 
the  static  stress-strain  curve.  C  is  the  same  as  before  and  H  is  deter¬ 
mined  by  higher  order  elastic  moduli.  Some  of  the  possibilities  are 
shown  in  Figure  6.  Either  positive  or  negative  solitary  waves  are 
possible  according  to  (6c) ,  and  (6d)  is  a  periodic  solution  of  unusual 
shape . 


w 


The  dynamic  trajectories,  D,  of  the  nonlinear  steady  waves  may 
deviate  considerably  from  the  static  curve  P(w',u,0)  =  0  in  the  w'-u 
plane.  To  see  this,  again  suppose  that  the  strain  energy  is  partitioned 
as  in  (42)  and  consider  the  case  A  ■  0.  Figure  (7a)  shows  a  contour 
plot  of  Wj  on  the  w'-u  plane  with  the  curve  P  =  Yl^  =0  and  several  lines 

of  u  =  const  drawn  in.  In  Figure  (7b)  several  curves  of  S  *  K.  vs. 

w' 

with  u  =  constant  are  shown.  Shown  as  a  dashed  curve  on  (7b)  is  the 
static  stress-strain  curve,  which  has  somewhat  smaller  slope  than  any  of 
the  curves  for  u  =  const,  according  to  (17).  From  reference  to  (58) .  it 
can  be  observed  that  points  of  the  nonlinear  trajectories  lie  on  the1 
intersection  of  the  surfaces  =  S(w',u)  and  S2  =  p^w'.  In  the 
projection  of  Figure  (7b)  is  a  straight  line  through  the  origin  with 
slope  PjC  ,  for  example  the^line  p^  (for  c  less  than  c^,  the  bar  speed 
at  w'  =  0).  By  comparing  the  line  PjP2  with  the  curves  for  S  at  u  *  const 
and  the  static  stress/strain  curve,  it  can  be  seen  that  between  points 
Pj  and  pj  the  dynamic  trajectory  deviates  to  values  of  u  less  than  those 
on  P  =  0,  and  at  points  outside  the  segment  p.p2  the  trajectory  deviates 
to  greater  values  of  u  than  those  on  P  =  0.  Thus,  the  trajectory,  D,  is 
qualitatively  like  the  dashed  curve  in  Figure  (7c) . 


Maximum  or  minimum  amplitudes  of  u  will  occur  where  u'  =  0.  With 
u'  =  0,  equation  (58)2  generates  curves,  C,  in  the  w'-u  plane  that 
correspond  to  the  intersection  of  the  surfaces  Wj  =  W^(w',u)  and 
^2  =  J  PjC^w'  -B.  The  intersection  of  curves  C  and  D  mark  the  ends  of 
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Figure  7:  Comparison  of  Steady  Dynamic  Trajectory  with  Static  Trajectory. 

a)  Contours  of  =  const,  and  Static  Curve  P=0.  b)  Intersec¬ 
tion  of  Surfaces  S(w',u)  and  p.c^w' .  c)  Dynamic  Trajectory 
Deduced  from  4b).  d)  Intersection  of  Cylinder  l/2p.c^w'^-B 
with  Widest  Projection  of  W,(w',u). 


the  dynamic  trajectory.  By  differentiating  (58)2  with  respect  to  w',  it 
may  be  seen  that  the  points  on  C  satisfy 

p  5^-  ♦  s  -  PjC2*'  .  (69) 

At  the  intersection  of  C  and  D,  (58)  ^  holds  as  well  so  that  at  those 
points 

P  3Vr  ■  o  .  (70) 

That  is,  either  P  =  0  or  the  slope  of  C  is  horizontal. 

Now  consider  further  the  intersection  of  surfaces  and  #2  with 
B  <  o  (B  >  o) .  Fig  (7d)  shows  the  projection  parallel  to  the  u  axis  of 
the  widest  points  of  =  W^(w',u).  These  points  correspond  to  the 
curve  P  ■  0  in  (7a),  and,  in  fact,  the  projected  potential  well  in  (7d) 

-  1  2  2 

generates  the  static  stress/strain  curve.  The  surface  “  2"  plc  W*  “  ® 

is  a  parabolic  cylinder  with  generators  parallel  to  the  u  axis.  An  end 
view  of  the  cylinder  is  shown  as  the  dashed  parabola  in  (7d) .  Since  the 
level  curves  of  are  closed,  as  shown  in  (7a),  and  since  the  level 
curves  of  W2  are  straight  lines  parallel  to  the  u-axis,  the  curves  C  are 
qualitatively  similar  to  those  sketched  in  (7c).  Thus,  for  the  case 
considered,  which  corresponds  to  Figure  (5d) ,  the  motion  lies  along  D 
and  either  oscillates  between  points  p^  and  p^  or  begins  at  Pj  and  moves 
to  the  right  until  interrupted  by  an  unsteady  motion.  In  either  case, 
note  the  divergence  of  D  from  the  curve  P  *  0. 

Other  cases  for  different  shapes  of  the  potential  well  W1(w,,u)  may 
be  considered  in  a  similar  manner.  However,  from  teh  preceding  qualita¬ 
tive  discussion  it  is  clear  that  te  cubic  expansion  in  equation  (63)  leads 
to  the  qualitatively  correct  curves  in  the  phase  plane  as  shown  in 
Figures  (4)  and  (5) ,  provided  that  the  function 

K(u'2)  =  Q(u')u'  -  W2(u')  -  \  P2c2  u'2 

2 

is  monotonic  and  hence  single  valued  in  u*  .  If  this  is  not  the  case, 
then  the  curves  of  uf^  vs.  u  in  Figures  (4),  (5),  and  (6)  will  have 
bifurcation  points.  The  analysis  would  then  be  somewhat  more  complicated 
and  will  not  be  attempted  here. 


VII.  REMARKS 

This  report  has  been  motivated  by  the  need  to  interpret  wave  propa¬ 
gation  experiments  in  rods  and  bars.  These  experiments  are  of  two  types. 
In  one,  a  cylindrical  bar  is  fired  from  an  air  gun  so  as  to  impact  an 
identical  bar,  which  is  at  rest  in  an  evacuated  chamber.  The  bar  at  rest 
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is  instrumented  at  several  stations  along  its  length  with  foil  resistance 
strain  gages.  In  the  other  experiment,  a  long  rod  penetrator,  which  is 
also  instrumented  with  strain  gages,  is  impacted  by  a  target  plate.  3 
In  either  case  a  transient  plastic  wave  propagates  into  the  bar,  strain¬ 
time  data  are  recorded  for  the  instrumented  stations  of  the  bar,  and 
then  the  data  are  interpreted  with  the  aid  of  equations  (2) .  By  compar¬ 
ing  with  streak  camera  records  Hauver  has  shown  that  (2)^  gives  very 
accurate  results  for  particle  velocity  ,  but  (2)2  always  gives  higher 
stresses  than  the  quasi-static  stress-strain  curve  even  though  the  high 
strength  steels  and  aluminums  used  are  not  rate  sensitive  in  Hopkinson 
bar  tests.  The  only  other  source  of  excess  stress  that  has  been 
suggested  is  the  effect  of  finite  transverse  dimensions,  hence  the 
present  interest  in  a  nonlinear  theory  that  accounts  for  radial  motion. 

This  report  has  dealt  only  with  nonlinear  elastic  materials,  not 
elastic-plastic  materials.  Nevertheless,  there  are  several  results  that 
are  at  least  suggestive  of  general  features  that  should  apply  to  the 
elastic-plastic  case.  These  are  as  follows. 

a.  The  averaging  in  equation  (10)  and  the  identifications 
made  in  equations  (11)  and  (12)  are  independent  of  constitutive  rela¬ 
tions.  Therefore,  the  equations  of  motion  (5)  hold  for  the  elastic- 
plastic  case  as  well  as  the  elastic  case,  but  of  course  equations  (6) 
will  not  be  appropriate  then. 

b.  With  v  ■  *  and  e  ■  w'  equations  (1)  are  entirely  consistent 
with  (5) .  As  a  consequence,  if  each  level  of  strain  propagates  at  a 
constant  velocity,  experimental  data  may  still  be  interpreted  through 
the  use  of  (2)  for  both  particle  velocity  and  stress.  In  this  integra¬ 
tion  no  use  is  made  of  (5)2  so  that  radial  strain  u  remains  undetermined. 
In  particular,  there  is  no  requirement  for  u  to  have  equilibrium  values 
corresponding  to  the  strain  e.  Thus,  the  computed  stress  S  should  be  as 
accurate  as  the  computed  particle  velocity  w  even  though  S  cannot  be 
measured  independently  as  can  w.  In  a  compressive  wave  the  excess 
stress  may  be  accounted  for  by  underexpansion  in  the  radial  direction 
caused  by  inertial  lag. 

c.  In  impact  experiments  of  the  type  described  here  it  is 
commonly  observed  that  the  plastic  wave  speed  becomes  nearly  constant 
at  the  higher  values  of  strain.  Thus,  steady  waves  are  of  interest  at 
least  as  asymptotic  solutions.  Because  of  dissipation  it  is  probably 
unreasonable  to  expect  that  solitary  waves  or  periodic  waves  will  occur 
in  elastic -plastic  materials,  but  pure  loading  waves  should  be  similar 
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in  either  elastic  or  elastic-plastic  materials.  Thus,  unbounded  waves 
of  the  type  shown  in  Figures  (4g),  (4h) ,  (5c),  or  (5d)  could  occur.  In 
an  actual  physical  problem  such  a  wave  must  be  interrupted  by  an 
unsteady  wave  or  perhaps  a  free  boundary  (as  in  the  penetration  experi¬ 
ment)  . 

Future  work  will  attempt  to  elaborate  on  some  of  these  issues. 
Initial -boundary  value  problems  will  be  posed  with  an  emphasis  on 
asymptotic  solutions.  The  elastic-plastic  case  will  also  be  formulated 
and  the  characteristic  features  of  those  solutions  examined  and  compared 
with  elastic  solutions. 


33 


REFERENCES 


1.  J.  F.  Bell,  The  Experimental  Foundations  of  Solid  Mechanics, 

Handbuch  der  Physilc,  Vol.  Vla/l,  Springer-Verlag,  New  York,  1973. 

2.  G.  E.  Hauver,  "Penetration  with  Instrumented  Rods,"  Int.  J.  Eng.  Sci., 
16,  1978,  pp.  871-877. 

3.  R.  Skalak,  "Longitudinal  Impact  of  a  Semi-Infinite  Circular  Elastic 
Bar,"  J.  Appl.  Mech. ,  24,  1957,  pp.  59-64. 

4.  W.  A.  Green,  "Dispersion  Relations  for  Elastic  Waves  in  Bars,"  in 
Progress  in  Solid  Mechanics,  Vol.  I,  ed.  I.  N.  Sneddon  and  R.  Hill, 
North-Holi and,  Amsterdam,  1960. 

5.  J.  H.  Shea,  "Propagation  of  Plastic  Strain  Pulses  in  Cylindrical 
Lead  Bars,"  J.  Appl.  Phys.,  39,  1968,  pp.  4004-4011. 

6.  G.  P.  DeVault,  "The  Effect  of  Lateral  Inertia  on  the  Propagation  of 
Plastic  Strain  in  a  Cylindrical  Rod,"  J.  Mech.  Phys.  Sol.,  13,  1965, 
pp.  55-68. 

7.  S.  S.  Ant man.  The  Theory  of  Rods,  Handbuch  der  Physik,  Vol.  VIa/2, 
Springer-Verlag,  New  York  1972. 

8.  J.  W.  Nunziato  and  E.  K.  Walsh,  "On  the  Influence  of  Void  Compac¬ 
tion  and  Material  Non-Uniformity  on  the  Propagation  of  One-Dimensional 
Acceleration  Waves  in  Granular  Materials,"  Arch .  Rat .  Mech .  Anal . , 

64,  1977,  pp.  299-316  and  Addendum,  Arch.  Rat.  Mech.  Anal.,  67,  1977, 
pp.  395-397. 

9.  J.  W.  Nunziato  and  E.  K.  Walsh,  "One-Dimensional  Shock  Waves  in 
Uniformly  Distributed  Granular  Materials,"  Int.  J.  Solids  and 
Structures,  14,  1978,  pp.  681-689. 

10.  D.  F.  Parker  and  B.  R.  Seymour,  "Finite  Amplitude  One-Dimensional 
Pulses  in  an  Inhomogeneous  Granular  Material,"  Arch.  Rat.  Mech.  Anal., 
72,  1980,  pp.  265-284. 

11.  M.  F.  McCarthy,  private  communication. 

12.  S.  S.  Antman,  "Qualitative  Theory  of  the  Ordinary  Differential 
Equations  of  Nonlinear  Elasticity,"  in  Mechanics  Today,  Vol.  1,  ed. 

S.  Nemat -Nasser,  Pergamon,  New  York,  1972. 

13.  R.  D.  Mindlin  and  G.  Herrmann,  "A  One-Dimensional  Theory  of  Com- 
pressional  Waves  in  an  Elastic  Rod,"  Proceedings  of  the  First  U.S. 
National  Congress  of  Applied  Mechanics,  1950,  pp.  187-191. 


34 


14.  G.  B.  Whitham,  "Some  Comments  on  Nave  Propagation  and  Shock  Nave 

Structure  with  Applications  to  Magnetohydrodynamics,"  Comm,  on  Pure 
and  Appl.  Math. ,  12,  1959,  pp.  113-158.  ~ 

15.  G.  B.  Whitham,  Linear  and  Nonlinear  Naves,  John  Niley,  New  York,  1974. 

16.  T.  T.  Nu,  "A  Note  on  the  Stability  Condition  for  Certain  Nave  Propa¬ 
gation  Problems."  Comm,  on  Pure  and  Appl.  Math.,  14,  1961,  pp.  745- 
747. 

17.  F.  Bowman,  "Introduction  to  Bessel  Functions,"  Dover,  New  York,  1958. 

18.  M.  Abromowitz  and  I.  A.  Stegun,  Handbook  of  Mathematical  Functions, 
Dover,  New  York,  1965. 

19.  J.  N.  Nunziato  and  E.  K.  Nalsh,  "Small -Amplitude  Nave  Behavior  in 
One-Dimensional  Granular  Solids,"  J.  Appl.  Mech.,  44,  1977,  pp.  559- 
564. 

20.  E.  Varley  and  E.  Cumberbatch,  "Non-Linear  High  Frequency  Sound  Naves," 
J.  Inst.  Maths.  Applies.,  2,  1966,  pp.  133-143. 

21.  S.  Leibovich  and  A.  R.  Seebass,  "Examples  of  Dissipative  and  Disper¬ 
sive  Systems  Leading  to  the  Burgers  and  the  Korteweg-de  Vries  Equa¬ 
tions,"  in  Nonlinear  Naves,  S.  Leibovich  and  A.  R.  Seebass,  eds., 
Corness  University  Press,  Ithaca,  N.Y.,  1974. 

22.  G.  A.  Nariboli  and  A.  Sedov, "Burgers's -Korteweg-de  Vries  Equation 
for  Visscoelastic  Rods  and  Plates,"  J.  Math.  Anal,  and  Appl.,  32, 

1970,  pp.  661-667. 

23.  G.  E.  Hauver  and  A.  Melani,  "Strain-Gage  Techniques  for  Studies  of 
Projectile  Behavior  During  Penetration,"  ARBRL-MR-03082,  Ballistic 
Research  Laboratory,  February  1981. 

24.  S.  C.  Chou,  AMMRC,  private  communication. 


35 


NEXT 

PAGE 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 


No.  of 

Copies  Organization 


1 2  Commander 

Defense  Technical  Info  Center 
ATTN:  DDC-DDA 
Cameron  Station 
Alexandria,  VA  22314 

1  Deputy  Assistant  Secretary 
of  the  Army  (RSD) 

Department  of  the  Army 
Washington,  DC  20310 

1  HQDA  (DAMA-ARP-P,  Dr.  Watson) 
Washington,  DC  20310 

1  HQDA  (DAMA-MS) 

Washington,  DC  20310 

1  Commander 

US  Army  Ballistic  Missile 
Defense  Systems  Command 
ATTN:  SENSC,  Mr.  Davidson 
P.  0.  Box  1500 
Huntsville,  AL  35804 

1  Director 

US  Army  Ballistic  Missile 
Defense  Systems  Office 
1320  Wilson  Boulevard 
Arlington,  VA  22209 

1  Director 

US  Army  Advanced  BMD 
Technology  Center 
ATTN:  CRDABH-5,  W.  Loomis 
P.0.  Box  1500,  West  Station 
Huntsville,  AL  35807 

1  Commander 

US  Army  War  College 
ATTN:  Lib 

Carlisle  Barracks,  PA  17013 


1  Commander 

US  Army  Command  and  General 
Staff  College 
ATTN :  Archives 
Fort  Leavenworth,  KS  66027 

1  Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCDMD-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

4  Commander 

US  Army  Armament  Research 
and  Development  Command 
ATTN:  DRDAR-TSS  (2  cys) 

DRDAR-SC ,  Dr.  E.  Bloore 
DRDAR-LC,  Dr.  J.  Frasier 
Dover,  NJ  07801 

l  Commander 

US  Army  Armament  Materiel 
Readiness  Command 
ATTN:  DRSAR-LEP-L,  Tech  Lib 
Rock  Island,  II  61299 

3  Director 

US  Army  Armament  Research 
and  Development  Command 
Benet  Weapons  Laboratory 
ATTN:  DRDAR-LCB-TL 

Dr.  F.  Schneider 
Dr .  T .  Davidson 
Watervliet,  NY  12189 

1  Commander 

US  Army  Aviation  Research 
and  Development  Command 
ATTN:  DRSAV-E 
P.  0.  Box  209 
St.  Louis,  MO  63166 


37 


distribution  list 


No.  of 

Copies  Organization 


No.  of 
Copies 


Organization 


Director 

US  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

Commander 

US  Army  Communications  Rsch 
and  Development  Command 
ATTN :  DRDCO-PPA-SA 
Fort  Monmouth,  NJ  07703 

l  Commander 

US  Army  Electronics  Research 
and  Development  Command 
Technical  Support  Activity 
ATTN:  DELSD-L. 

Fort  Monmouth,  NJ  07703 

l  Commander 

US  Army  Harry  Diamond  Labs 
ATTN:  DELHD-TA-L 
2800  Powder  Mill  Road 
Adelphi,  MD  20783 

1  Commander 

US  Army  Missile  Command 

ATTN:  DRSMI-R 

Redstone  Arsenal,  AL  35809 

1  Commander 

US  Army  Missile  Command 
ATTN:  DRSMI-YDL 
Redstone  Arsenal,  AL  35809 


L  Commander 

US  Army  Napick  Research 
and  Development  Center 
ATTN:  DRXRE,  Dr.  D.  Sieling 
Natick,  MA  01762 

1  Commander 

US  Army  Tank  Automotive 

Research  8  Development  Cmd 
ATTN :  DRDTA-UL 
Warren,  MI  48090 

1  Commander 

US  Army  Electronics  Proving 
Ground 

ATTN :  Tech  Lib 

Fort  Huachuca,  AZ  85613 

3  Commander 

US  Army  Materials  and 

Mechanics  Research  Center 
ATTN:  DRXMR-T,  J.  Mescal 1 

R.  Shea 

S.  C.  Chou 

Watertown,  MA  02172 

3  Commander 

US  Army  Research  Office 
ATTN:  Dr.  G.  Mayer 

Dr.  F.  Schmiedeshoff 
Dr.  E.  Saibel 
p.  0.  Box  12211 
Research  Triangle  Park 
NC  27709 


2  Commander 

US  Army  Mobility  Equipment 
Research  8  Development  Cmd 
ATTN:  DRDME-WC 
DRSME-RZT 

Fort  Belvoir,  VA  22060 


Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL,  Tech  Lib 
White  Sands  Missile  Range 
NM  88002 


38 


DISTRIBUTION  LIST 


No .of  No .  of 

Copies  Organ! zat ion  Copies  Organization 

1  Office  of  Naval  Research  1  RADC  (EMTLD,  Lib) 

Department  of  the  Navy  Griff iss  AFB,  NY  13440 

ATTN:  Code  402 

Washington,  DC  20360  1  AUL  (3T-AUL-60-118) 

Maxwell  AFB,  AL  36112 

1  Commander 

Naval  Surface  Weapons  Center  1  AFFDL/FB  (Dr.  J.  Halpin) 

ATTN:  Code  GR-9,  Dr.  W. Soper  Wright-Patterson  AFB,  OH  45433 

Dahlgren,  VA  22448 

1  AFML  (Dr.  T.  Nicholas) 

1  Commander  Wright-Patterson  AFB,  OH  45433 

Naval  Weapons  Center 


ATTN:  Code  3813,  M.Backman 

1 

Director 

China  Lake,  CA  93557 

Lawrence  Livermore  Laboratory 
ATTN:  Dr.  M.  Wilkins 

Commander  and  Director 

Livermore,  CA  94550 

Naval  Electronics  Laboratory 

ATTN:  Lib 

4 

Sandia  National  Laboratories 

San  Diego,  CA  92152 

ATTN:  Dr.  L.  Davison 

Dr.  P.  Chen 

Commander 

Dr.  L.  Bertholf 

Naval  Research  Laboratory 

Dr.  W.  Herrmann 

ATTN:  Code  5270,  F. MacDonald 

Albuquerque,  NM  87115 

Code  2020,  Tech  Lib 

Code  7786,  J.  Baker 

1 

IBM  Watson  Research  Center 

C.  Sanday 

ATTN:  R.A.  Toupin 

Washington,  DC  20375  Poughkeepsie,  NY  12601 


1  AFATL  (DLDG) 

Eglin  AFB,  FL  32542 

1  AFATL  (DLDL,  MAJ  J.  D.  Morgan) 
Eglin  AFB,  FL  32542 

1  AFATL  (DLYW) 

Eglin  AFB,  FL  32542 

1  AFATL  (J.  Smith) 

Eglin  AFB,  FL  32542 


5  Brown  University 

Division  of  Engineering 
ATTN:  Prof.  R.  Clifton 
Prof.  H.  Kolsky 
Prof.  A.  Pipkin 
Prof.  P.  Symonds 
Prof.  J.  Martin 
Providence,  RI  02192 

3  California  Institute  of  Tech 
Division  of  Engineering  and 
Applied  Science 
ATTN:  Dr.  J.  Miklowitz 
Dr.  E.  Sternberg 
Dr.  J.  Knowles 
Pasadena,  CA  91102 


39 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

4  Carnegie  Mellon  University 
Department  of  Mathematics 
ATTN:  Dr.  D.  Owen 

Dr.  M.  E.  Gurtin 
Dr.  B.  Coleman 
Dr.  W.  Williams 
Pittsburgh,  PA  15213 

2  Catholic  University  of  America 
School  of  Engineering  and 
Architecture 
ATTN:  Prof.  A.  Durelli 
Prof.  J.  McCoy 
Washington,  DC  20017 

1  Harvard  University 
Division  of  Engineering 

and  Applied  Physics 
ATTN:  Dr.  G  Carrier 
Cambridge,  MA  02138 

2  Iowa  State  University 
Engineering  Research  Lab 
ATTN:  Dr.  G.  Nariboli 

Dr.  A.  Sedov 
Ames,  IA  50010 

4  Johns  Hopkins  University 
ATTN:  Prof.  J.  Bell 

Prof.  R.  E.  Green 
Prof.  C.  A.  Truesdell 
Prof.  J.  L.  Ericksen 
34th  and  Charles  Streets 
Baltimore,  MD  21218 

2  Lehigh  University 

Center  for  the  Application 
of  Mathematics 
ATTN:  Dr.  E.  Varley 
Dr.  R.  Rivlin 
Bethlehem,  PA  18015 


No.  of 

Copies  Organization 

1  Massachusetts  Institute  of 
Technology 

ATTN:  Dr.  R.  Probstein 
77  Massachusetts  Avenue 
Cambridge,  MA  02139 

1  New  York  University 
Department  of  Mathematics 
ATTN:  Dr.  J.  Keller 
University  Heights 
New  York,  NY  10053 

2  North  Carolina  State  University 
Department  of  Engineering 

Mechanics 

ATTN:  Dr.  W.  Bingham 
Dr.  Y.  Horie 
P.  0.  Box  5071 
Raleigh,  NC  27607 

1  Pennsylvania  State  University 
Engineering  Mechanical  Dept. 
ATTN:  Prof.  N.  Davids 
University  Park,  PA  16502 


1  Purdue  University 
Institute  for  Mathematical 

Sciences 

ATTN:  Dr.  E.  Cumberbatch 
Lafayette,  IN  47907 

2  Rice  University 
ATTN:  Dr.  R.  Bowen 

Dr.  C.  C.  Wang 
P,  0.  Box  1892 
Houston,  TX  77001 


40 


DISTRIBUTION  LIST 


No.  of 
Copies 


Organization 


No.  of 
Copies 


Organization 


1  Southern  Methodist  University  1 
Solid  Mechanics  Division 

ATTN:  Prof.  H.  Watson 
Dallas,  TX  75221 

2  Southwest  Research  Institute 

Department  of  Mechanical  3 

Sciences 

ATTN:  Dr.  U.  Lindholm 
Dr.  W.  Baker 
8500  Culebra  Road 
San  Antonio,  TX  78228 

1  Temple  University 

College  of  Engineering  Tech.  2 

ATTN:  Dr.  R.  M.  Haythomthwaite, 

Dean 

Philadelphia,  PA  19122 

1  Tulane  University 

Dept  of  Mechanical  Engineering 
ATTN:  Dr.  S.  Cowin  1 

New  Orleans,  LA  70112 

2  University  of  California 
ATTN:  Dr.  M.  Carroll 

Dr.  W.  Goldsmith 

Berkeley,  CA  94704  2 

1  University  of  California 
Dept,  of  Aerospace  and 
Mechanical  Engineering  Science 
ATTN:  Dr.  Y.  C.  Fung 
P.  0.  Box  109 
La  Jolla,  CA  92037 


University  of  Delaware 
Department  of  Mechanical 
Engineering 

ATTN:  Prof.  J.  Vinson 
Newark,  DE  19711 

University  of  Florida 
Dept,  of  Engineering  Science 
and  Mechanics 

ATTN:  Dr.  C.  A.  Sciammarilla 
Dr.  L.  Malvern 
Dr.  E.  Walsh 
Gainesville,  FL  32601 

University  of  Houston 
Department  of  Mechanical 
Engineering 
ATTN:  Dr.  T.  Wheeler 

Dr.  R.  Nachlinger 
Houston,  TX  77004 

University  of  Illinois 
Dept,  of  Theoretical  anc 
Applied  Mechanics 
ATTN:  Dr.  D.  Carlson 
Urbana,  IL  61801 

University  of  Illinois 
at  Chicago  Circle 
College  of  Engineering 
Dept,  of  Materials  Engineering 
ATTN:  Prof.  A.  Schultz 
Dr.  T.  L.  T.  Ting 
P.  0.  Box  4348 
Chicago,  IL  60680 


1  University  of  California  1  University  of  Iowa 

Department  of  Mechanics  ATTN:  Dr.  L.  Valanis 

ATTN:  Dr.  R.  Stem  Iowa  City,  IA  52240 

504  Hilgard  Avenue 
Los  Angeles,  CA  90024 


41 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

4  University  of  Kentucky 
Department  of  Engineering 
Mechanics 

ATTN:  Dr.  M.  Beatty 

Prof.  0.  Dillon,  Jr. 
Prof.  P.  Gillis 
Dr.  D.  Leigh 
Lexington,  KY  40506 

1  University  of  Maryland 
Department  of  Mathematics 
ATTN:  Dr.  S.  Antman 
College  Park,  MD  20742 

2  University  of  Maryland 
Department  of  Mechanical 

Engineering 
ATTN:  Prof.  J.  Yang 
Dr.  J.  Dally 
College  Park,  MD  20742 

1  University  of  Minnesota 
Department  of  Engineering 
Mechanics 

ATTN:  Dr.  R.  Fosdick 
Minneapolis,  MN  55455 

1  University  of  Notre  Dame 
Department  of  Metallurgical 
Engineering  and  Materials 
Sciences 

ATTN:  Dr.  N.  Fiore 
Notre  Dame,  IN  46556 

1  University  of  Pennsylvania 
Towne  School  of  Civil  and 
Mechanical  Engineering 
ATTN:  Prof.  Z.  Hashin 
Philadelphia,  PA  1910S 


No.  of 

Copies  Organization 

4  University  of  Texas 

Department  of  Engineering 
Mechanics 

ATTN:  Dr.  M.  Stern 

Dr.  M.  Bedford 
Prof.  Ripperger 
Dr.  J.  T.  Oden 
Austin,  TX  78712 

1  University  of  Washington 
Department  of  Mechanical 

Engineering 

ATTN:  Prof.  J.  Chalupnik 
Seattle,  WA  98105 

2  Yale  University 
ATTN:  Dr.  B.  Chu 

Dr.  E.  Onat 
400  Temple  Street 
New  Haven,  CT  96520 

2  Washington  State  University 
Department  of  Physics 
ATTN:  Dr.  R.  Fowles 
Dr.  G.  Duvall 
Pullman,  WA  99163 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 

ATTN:  DRXSY-MP,  H.  Cohen 
Cdr,  USATECOM 

ATTN:  DRSTE-TO-F 
Dir,  USACSL ,  Bldg.  E3S16,  EA 
ATTN:  DRDAR-CLB-PA 


42 


USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3,  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc . )  _ _ _ _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/ contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 


Name: 

Telephone  Number: 
Organization  Address: 


—  -  -  -  -  -  -  -  - FOLD  HERE  - 


Director 

US  Army  Ballistic  Research  Laboratory 
Aberdeen  Proving  Ground,  MD  21005 


NO  POSTAGE 
NECESSARY 
IF  MAILED 
IN  THE 

UNITED  STATES 


OFFICIAL  BUSINESS 

PENALTY  FOR  PRIVATE  USE,  *300 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  NO  12062  WASHINGTON.OC 
POSTAGE  WILL  BE  PAIO  BY  DEPARTMENT  OF  THE  ARMY 


Director 

US  Army  Ballistic  Research  Laboratory 
ATTN:  DRDAR-TSB 

Aberdeen  Proving  Ground,  MD  21005 


FOLD  HERE 


